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PhI Abstract 

We present a comprehensive analysis of the Non-Minimal Supersymmetric Stan- 



dard Model (NMSSM) for large values of tan f3, the ratio of the vacuum expectation 



^ ■ values of the two Higgs doublets, which arise when we impose the constraint of 

the unification of Yukawa couplings in the model. In this limit we show that the 
vacuum expectation value of the singlet is forced to be large, of the order of 10 TeV. 
The singlet decouples from the lightest CP-even neutral Higgs boson and the neu- 
tralinos. We compare our results with the corresponding particle spectrum of the 
Minimal Supersymmetric Standard Model in the same limit. With the exception of 
the lightest Higgs boson, the particle spectrum in the model turns out to be heavy. 
The Higgs boson mass, after the inclusion of radiative corrections, is found to be in 
the neighbourhood of ~ 130 GeV. 
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1 Introduction. 



The recent determination of the couphng constants of the SU{3)c x SU{2)l x U{1)y 
model of strong and electroweak interactions is compatible with the super symmetric uni- 
fication of these couplings with a supersymmetry 0] breaking scale of the order of 1 
TeV. In supersymmetric (SUSY) models, at least two higgs doublets {Hi and H2) with 
opposite hypercharge are required to give masses to the up- and down- type quarks (and 
charged leptons), and to cancel gauge anomalies. The ratio of the vacuum expection values 
of the neutral components (< >= vi and < H2 >= V2) of these two Higgs doublets, 
tan/3 = ^, is a crucial parameter for the predictions from supersymmetric grand uni- 
fied theories. One particularly predictive framework is based on the assumption that the 
heaviest generation fermions lie in a unique 16- dimensional representation of the unifying 
gauge group 5*0(10), with the Higgs doublets in a 10- dimensional representation of the 
group [3]. Furthermore, if one makes the additional assumption that the fermion masses 
are generated by a single complex 10-plet of SO{10) through a /i. 16. 16. 10 term in the 
superpotential at the unification scale Mx, determined from gauge coupling unification, 
one would have the prediction 0: 

ht = hb = hr = h (1-1) 

where ht^b,T are the Yukawa couplings of t, b and r. The coupled system of gauge and 
Yukawa couplings are then evolved from Mx down to the weak scale, and tan jS is deter- 
mined from the accurately measured value of rrir = 1.78 GeV. When h is chosen in such 
a manner as to yield a value for mi,{mb) in its observed range of 4.25 ± 0.10 GeV 0, a 
rather good prediction for the top-quark mass is obtained, which, with the present central 
value of as{Mz) = 0.12, lies in the range favored by the experimental data In such a 
situation, tan P is found to saturate what is considered to be a theoretical upper bound 
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on its values of mt/mf,, and the Yukawa coupling h is found to come out to be rather large 
0(1 — 3), with a certain insensitivity to the exact value, since it is near a fixed point of 
its evolution. 

At the weak scale, the minimal particle content of such a grand unified theory is 
that of the minimal supersymmetric standard model with each bosonic (fermionic) 
degree of freedom of the standard model complemented by fermionic (bosonic) one, apart 
from having two Higgs doublets as mentioned earlier. However, it does not involve in 
any great detail the remaining aspects of the embedding of the standard model into a 
supersymmetric grand unified framework. With the particle content of the MSSM and 
an additional discrete symmetry, called R-parity which forbids couplings that can 
lead to rapid nucleon decay, it is possible to construct a self-consistent and successful 
framework. In particular, the Higgs sector of the MSSM contains five physical degrees of 
freedom: two CP even {h^ and H^) and one CP-odd (A°) neutral and complex charged 
Higgs bosons (i?^). The spectrum of the Higgs bosons is strongly affected by the grand 
unified theory assumptions such as the ones described above. For instance, the mass 
of the lightest higgs boson m/jo, which is related through the D-term in the potential 
to the mass of the Z-boson (m^), after inclusion of radiative corrections is found to 
be ~ 140 GeV 0. Furthermore, supersymmetry breaking is understood to arise from 
embedding the MSSM into a supergravity framework and writing down all the possible 
soft supersymmetry breaking terms consistent with gauge and discrete symmetries that 
define the model. It is ususally assumed that many of the parameters describing these 
terms are in fact equal at the unification scale in order to have a predictive framework, 
which is motivated by such arguments as the weak principle of equivalence when applied 
to coupling of the SUSY breaking hidden sector to the known sector via gravitation. 
Such universal boundary conditions have had to be modified somewhat in order to avoid 
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fine tuning and to minimize the effects of possible radiative corrections to the b-quark 
mass induced through the supersymmetry breaking sector of the theory at the one-loop 
level 

In SU{5) type unification where tan (3 is free, the region tan/? ~ 1 is also a region 
which is favoured for the unification of the b-quark and r-lepton Yukawa couplings from 



the observed data [T^. One crucial difference between the two extremes is that for large 
values of tan/?, the Yukawa couplings of the b-quark (and that of the r-lepton) always 
remain comparable to that of the top-quark, with the observed hierachy in their masses 
arising from the large value of tan (3; when tan /5 ~ 1, the Yukawa couplings of the b-quark 
and that of the r-lepton are much smaller than that of the top-quark. 

Despite its many successes, it may be premature to confine our attention only to the 
MSSM, especially because of the presence of dimensionful Higgs bilinear parameter /z in 
the superpotential. An alternative to the MSSM that has been widely considered is the 
Non-Minimal Supersymmetric Standard Model (NMSSM) , where the particle content of 
the Minimal Supersymmetric Standard Model (MSSM) is extended by the addition of 
a gauge singlet chiral superfield S, and one in which dimensionful couplings are eliminated 
through the introduction of a discrete symmetry. While the explicit coupling fiHiH2 is 
forbidden, an effective /i term is generated by the vacuum expectation value < S' > (= s) 
of the singlet Higgs field. 

The NMSSM is characterised by the superpotential [|12 



W = htQ- H2t% + hQ ■ H^b% + KL ■ H^t'^ + XSH1.H2 + \kS\ (1.2) 

where we have written the interactions of only the heaviest generation and the Higgs 
sector of the theory. In addition, one must add to the potential obtained from eq.( p..2| ), 
the most general soft supersymmetry breaking terms. These are: 

{htAtQ ■ H2t\ + hQ.Hib% + KArL ■ Hirf, + XA^Hi ■ H2S + hAkS^) + h.c. 



(1.3) 

in the standard notation. In addition, there are soft SUSY breaking mass terms for U (l)y , 
SU{2)l and SU{?>)c gauginos (A'; A", a = 1, 2, 3; A*; i = 1, ••••,8), which we shall denote 
by Ml, M2 and M3, respectively. We note that if = 0, the Lagrangian obtained from 
(1.2) has a global U{1) symmetry corresponding to — A^e*^, H1H2 HiH2e~^^^ , which 
is broken by the vacuum expectation values (VEVs) of the Higgs fields. In order to avoid 
axion associated with this symmetry we require k ^ 0. 

In this paper we present a study of the particle spectrum of the NMSSM, characterized 
by eq.( p..2|) and eq.( |1.3|) , for large values of tan/3 in detail, and compare it with the 
corresponding spectrum in the MSSM. We carry out a renormalization group analysis of 
this model with universal boundary conditions and analyze the RG improved tree-level 
potential at the scale Qq. The cut-off scale for the renormalization group evolution is 
chosen to be the geometric mean of the scalar top quark masses which is roughly equal 
to the corresponding mean of the scalar b-quark masses as well, since the evolution of 
the Yukawa couplings of the t- and b-quarks are equal upto their hypercharges and the 
relatively minor contribution of the r-lepton Yukawa coupling [we note here that this 
feature is crucial in justifying our use of the tree-level potential since with large tan (3, 
the Yukawa coupling of the b-quark is large, influencing the contributions of the large 
logarithms in the one-loop potential]. Whereas in the MSSM the parameters /i and B (the 
soft SUSY breaking parameter corresponding to the bilinear term in the superpotential 
of the MSSM) do not enter into the evolution of the other parameters of the model at the 
one-loop level, the situation encountered here is drastically different with a systematic 
search in the parameter space having to be performed with all parameters coupled from 
the outset. Our analysis of the minimization conditions that ensure a vacuum give rise to 
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severe fine tuning problems that are worse in the NMSSM in comparison to those in the 
MSSM. The problem is further compounded by our having to satisy three minimization 
conditions in contrast to the two that occur in the MSSM. In studies of the model where 
tan /3 is a free and adjustable parameter, the tuning of parameters is possible in order to 
meet all the requisite criteria, viz., minimization conditions, requirement that the vacuum 
preserve electric charge and color, etc. However in the present case where tan (3 is fixed 
and large, what we find is a highly correlated system. 

The plan of the paper is as follows: In Section 2, we review the basic framework of 
the MSSM at large tan (3 so as to set up the stage for our analysis of the NMSSM in the 
same limit, and thereby enable us to compare and contrast the two models. In Section 3 
a detailed analysis of the NMSSM in this limit is presented. An important conclusion of 
our analysis is that the vacuum expectation value of the singlet is forced to be large in 
this limit. In Section 4 we present a numerical study of the model and the conclusions. 



2 MSSM at Large tan/3 

As mentioned in the introduction, in the minimal supersymmetric extension of the stan- 
dard model, the ratio of the vacuum expectation values of the two Higgs doublets, tan /3, 
is an important parameter that determines the spectrum in an essential way. The unifi- 
cation of the third family Yukawa couplings, eg . ( |1 . 1|) , leads to tan/5 being determined in 



the MSSM, thus, reducing the parameter space significantly. The analysis of the MSSM 
starts with the scalar potential, where the parameters are evolved according to the one- 
loop renormalization group equations from Mx to the low energy scale Qo • At tree- level 
it is possible to analytically find the minimum of the potential, and we simply quote the 
results here [|14|. We start with the definitions: 
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where fnfj^, i = 1,2 are the soft-supersymmetry breaking mass parameters for the Higgs 
fields, and B is the soft-supersymmetry breaking bihnear parameter corresponding to the 
fiHiH2 term in the superpotential of MSSM. In order that the minimum of the potential 
break SU{2) x f/(l) gauge symmetry to U{l)em, we must have fiffil < whereas to 
prevent the potential from being unbounded from below we require + ^\ > 2|/i3p (note 
that the left-hand side of this is nothing but the square of the mass of the CP-odd Higgs 
boson, m^). These requirements tend to favor a situation when only one of /if or is 
negative, or where both are positive and one is somewhat smaller than the other. From 
here we see that radiative electroweak symmetry breaking via the effects of ht may be 
implemented in the MSSM in general as it drives to smaller (possibly negative) 



values in relation to fif, while /Xg can be easily set negative. 

In the attractive scenario of radiative electroweak symmetry breaking, the mass pa- 
rameters nl and III start out at the GUT scale with a universal positive value 

2 2 2 2 

= ^^ = mQ + fx , 

where ttiq is the universal scalar mass. Thus the symmetry is not broken at this scale. 
However, in the RG evolution to the electroweak scale, the large Yukawa coupling of the 
top quark to H2, which gives the top its mass, also drives the mass-squared parameter 
fi2 of H2 negative (competing against the QCD coupling), while the absence of a large 
Yukawa in the down sector keeps the mass squared of Hi positive. The above conditions 
are easily satisfied for a large range of initial conditions if h ~ 0(1), resulting in a very 
natural picture of electroweak symmetry breaking. This picture is essentially lost in the 
large tan P scenario, for the reasons discussed below. 

Firstly, since all the Yukawa couplings are comparable, the two higgs doublets tend 
to run in the same way, so either both /if and /ij stay positive at the electroweak scale 
and the symmetry does not break, or both become negative and the potential becomes 



unbounded from below. The effects whicfi differentiate between tfieir running, sudi as 
differing hypercharge assignments, are small, and a poor replacement for the usual ht ^ 
condition. Furthermore, an 0(1) splitting between ht and hb^r here is still of little use 
since it is quickly reduced in importance by the proximity to the fixed point. 

Secondly, even when the electroweak symmetry is broken, a large hierarchy of VEVs 
must be generated between the two similarly evolving Higgs fields. To see the implications 
of this, let us recall that the minimization of the tree level potential yields 

H-l -|- 

Prom here we see that large tan/? will require ix\ to be small in magnitude relative to 
+ Ii2- The second minimization condition is 

tan2/j = 4^^, (2.2) 

or, equivalently. 



2 _ 2(/xj-/x|tan^ (3) _ 2{m\^ - m]j^tiiv? „ o\ 

tan2/3-l - tan2/3-l ' ^^'^^ 



which in turn implies that for tan/? 1 the solution we are looking for requires ~ 
— m|/2, thus setting the scale. We, thus, see that a large hierarchy of VEVs requires 
the large hierarchy A'-a -C A'-i + A*2- This, as is well known, implies a degree of fine tuning 
between some of the parameters in the Lagrangian. Knowledge of m^^^ and m^^^ Qo 
determines /x^((5o) from the approximate relation 

,^(Q„) . - f ^ (2.4) 

2 

and then B = ^ \% also determined. Note that |//|| has to be small at large values 
of tan/3, so that B must also be small. In numerical studies it has been found that 
electroweak symmetry works well for reasonably large ht so long as h-^ is neglected. Upon 



its inclusion, the realization of electroweak symmetry breaking engenders find tuning. The 
dependence of /i is somewhat complicated, although it may seem from eq. that small 
fi gives large tan (3. On the other hand A (the trilinear soft SUSY breaking parameter) 
seems to be relatively less important in absolute terms (i.e., one does not obtain any 
significant restriction on its magnitude), although it does infiuence the running of the 
soft-parameters. In the NMSSM on the other hand, it plays an important role as we 
shall see in the next section. We note here that the MSSM spectrum is unchanged if 
Ml/2 —Ml/2 (where M1/2 is the universal gaugino mass) when A = 0. 

The main features of the sparticle spectrum turn out to be governed by (a) strong 
fi — Ml/2 correlation, (b) large values of M1/2 ~ 400 GeV, and (c) M1/2 ~ mo (uni- 
versal soft SUSY breaking scalar mass). This implies small mixing in the chargino and 
neutralino sectors. The lightest supersymmetric particle is mainly a bino, with a mass 
~ 0.4 Ml/2. The second lightest neutralino and lightest chargino are winos and hence 
almost degenerate in mass. The heaviest neutralino and chargino are Dirac (pseudo-Dirac 
in the case of neutralino) particles, with masses approximately equal to the parameter 
\fi{Qo)\. An important tree level relation is obtained for large tan/5 > 10, for which the 
tree-level value of the lightest CP-even higgs mass {rrih) is equal to rriz, whenever the CP- 
odd boson mass (m^) is larger than mz, while for itla < friz, rrih = mA- This tree level 
relation is approximately stable under radiative corrections, with the only difference that 
the range for which rrih = f^A holds, extends to values of tua somewhat larger than mz- 
Therefore, large values of tan (3 and values of CP-odd Higgs in the desired range imply 
iTih < mz- We note that large values of tan/3 and values CP-odd higgs mass niA < 70 
GeV are preferred to improve the agreement with the value of Rb = r(z-^Tdrons) iiieasured 
at LEP [|1|]. 
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3 NMSSM at Large tan P 

The potential for the Higgs fields of the NMSSM can be obtained from eq.(|l]2|) and eq.(|L3| 



through a standard procedure |12|. The minimization conditions (evaluated at Qo after 
all the parameters are evolved via their one- loop RG equations down to this scale), that 
determine the soft SUSY breaking Higgs masses in terms of the other parameters, are 



= -A^s(A, + ks) - X\vl + s') + + g''){vl - v^), (3.1) 

= -A-s(A, + ks) - X\vf + s') + -W + g'^){vl - vj), (3.2) 

f 2 4 

^1 = _A2(t;2 + vl) - 2^3'' - 2Xsv,V2 - kAkS - (3.3) 



The first two minimization conditions can be rewritten as: 



tan (3 = 2 , x2 2 ' (3-4) 



-2Xs){Ax + ks) 
ml^ +ml^+X\2s^ + v^) 



sin 2/3 = —2 — —2 — , ,^,0.2 , ..2x - (3-5) 



[Our normalization is such that v'^ = vf + f|(= 174GeV)2 and m| = ^{g'^ + g''^)v^, where 
g and g' are the gauge couplings of SU{2) and f/(l), respectivly.] These two equations 
give us some insight into the manner in which our solutions are likely to behave. Eq. 
( p.4|) guarantees that tan/? must lie between 1 and mt/mb. The proof, as in the case 
of the MSSM, relies once more on the RG equations that govern the behaviour of mass 
parameters and may be proved simply by reductio ad absurdum. For this purpose we 
need only consider the following equation expressing the momentum dependence of the 
difference of two supersymmetry breaking scalar mass parameters: 

- <J = + ^^1^^ + h'r^r) (3.6) 
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where t = log(/i), the logarithm of the momentum scale, and Xj, i = t, b, t, are combi- 
nations of scalar masses and trilinear couplings: 

Xt = ml + ml + ml^ + A^t, 

Xb = ml + ml + ml^ + Al (3.7) 
X^ = ml + ml + m]j^ + A^. 

It must be noted that in order to prove that tan/3 > 1 we neglect hi, and h^, and for 
proving tan/? < mt/mj, we retain them. 

From eg. ( p.4|) it is clear that, as in MSSM, in order to have large tan/? with m|^^ 
and m'jj^ being essentially degenerate, because the top and bottom Yukawa couplings are 
comparable in the RG evolution, the denominator of the equation has to be small at the 
weak scale. This implies the fine tuning condition 

m^^ + AV f« -m|/2. (3.8) 

From this condition it follows that the correspondence with the MSSM will occur in a 
certain well defined manner with the identification of As with /i. Similarly, we must 
identify Ax + ks with B, the bilinear soft supersymmetry breaking parameter of the 
MSSM. We will show below that for large values of tan/5 this identification occurs in a 
novel manner, not generic to the model, say, for tan/3 ~ 1. For large values of tan/3, 
eq.(|3.5|) implies 

Ax ^ -ks (3.9) 

This is analogous to the condition B ^ of the MSSM. The situation here is complicated 

because Ax is not a parameter that is fixed at Qq but is present from the outset. This is 

the first of the fine tuning problems we encounter in the NMSSM. 

A rearrangement of eq. (|3.5| ) equation yields 

X / . , N ^/ 9 ,a o, ml (tan/3^ - 1) A^f^sin2/3 

Xs{Ax + ks) = tan/3(-m^^ - AV)^ ^ (3.10) 
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In this equation we can discard the last term when tan j3 is large. Then, with the identifi- 
cation of the appropriate parameters in terms of those of the MSSM as described earlier, 
we recover all the MSSM relations of the previous sections without having to go through a 
limiting procedure that is required in the general case. The third minimization condition, 
eq.(|3.3|), can be rewritten as: 



m| = -AV-2A2 + _iL^ + AAAfc + ^. (3.11) 

In order to satisfy this condition, one must have large cancellations between the fourth 
and the first two terms, since the terms proportional to sin 2(3 are negligible. This requires 
that Ax and come out with the same sign and that their product be sufficiently large. It 
is this fine tuning condition that leads to problems with finding solutions with sufficiently 
small trilinear couplings in magnitude [0. 

The starting point of our analysis of the NMSSM is the estimation of the scale Mx with 
the choice of the SUSY breaking scale Qo ~ 1 TeV . For asijnz) = 0.12, Qq = 1 TeV and 
a = 1/128, we find upon integrating the one- loop beta functions, Mx = 1.9 x 10^^ GeV, 
and the unified gauge coupling adMx) = 1/25.6. We then choose a value for the unified 
Yukawa coupling h of 0(1). The free parameters of the model are the common gaugino 
mass (M1/2), the common scalar mass (mo), the common trilinear scalar coupling (A), 
and the two additional Yukawa coupligs (A, k), respectively. Note that our convention 
requires us to choose A > and A; < in order to conserve CP in the Yukawa sector of 



the model |]T3|. We also impose the constraint that \A\ < 3mo in order to guarantee 
the absence of electric-charge breaking vacua. In the case at hand this choice may have 
to be strengthened further due to the presence of large Yukawa couplings for the b-quark. 
The situation is considerably less restrictive when mild non-universality is allowed and, 
for instance, if strict Yukawa unification is relaxed. Given these uncertainties, we choose 
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to work with this constraint. We also compute the mass of the charged Higgs boson [|T2 



ml = ml, - AV - \{A^ + fcs)-^, (3.12) 

sm 2p 

where m^^ is the mass of the W-boson. We note that the radiative corrections to the 
charged Higgs mass are small for most of the parameter range, as in the case of MSSM |l^ . 



The reason for this is that a global SU{2) x SU{2) symmetry [IS] protects the charged 
Higgs mass from obtaining large radiative corrections. If this quantity were to come out 
to be negative, then the resulting vacuum would break electric-charge spontaneously and 
the corresponding point in the parameter space would be excluded. 

Choosing a particular set of values for the input parameters (M1/2, fno, A, A, k), 
satisfying the above constraints, we calculate the crucial parameters tan/?, r(= s/v), Ax 
and Ak and ks, that determine the physical spectrum, from their renormalization group 
evolution to Qq. We choose these input parameters in such a manner so as to ensure 



that we are in the neighborhood of a vacuum that breaks SU{2) x f/(l) [|l^]. With 
the boundary condition (1.1) at the unification scale, we find IM1/2I ~ 500 GeV for a 
wide range of values of A ~ 0.4 — 0.6 with k = —0.1, leading to values of tan/3 ~ 60. 
These magnitudes emerge when we choose A in a manner that does not lead to very 
small values of k, and furthermore, we note that the solutions depend only on the ratio 
X/k [|13|. Once the sign of k is fixed [from requirements of CP invariance], a solution is 
found only when the sign of M1/2 is negative. Phenomenological requirements enforce 
mo also to be rather large in order to guarantee that the lighter scalar tau be heavier 
than the lightest neutralino. The minimization condition (3.3) or its equivalent (3.11) 
enforces a large ratio [A/mo[ ~ 3. It then follows from eq.(3.9) that s must be rather 
large in comparison with mz, the only other scale in the problem. We note that here we 
are near the fixed point of A, and as such we expect that this provides us with a lower 
bound on the singlet vacuum expectation value s. Essentially, this implies that in order 
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to have Yukawa unification (1.1) in the NMSSM, and the resulting large value of tan/3, 
the singlet vacuum expectation value must be large compared to the doublet vacuum 
expectation values. As we shall see in Sec. 4, these conclusions are borne out by the 
detailed numerical calculations. 

We now analyze the implications of the large values of the singlet vacuum expectation 
value s, that arises in the NMSSM with eq. (1.1), on the mass spectra. In the basis 
(Imi/?, ImH^, ImS), with neutral Goldstone boson G° = cos P{ImH^) - sin P{lmH^) and 
the orthogonal combination = cos /3(Imi7°) +sin /3(Imi?2)) the tree-level mass squared 
matrix for the two pseudoscalar neutral Higgs fields can be written as (in numerical 
calculations we take into account one-loop radiative corrections in the effective potential 
approximation) 



^ "Mf XviA^ + 3ks) ^ 



y Xv{Aj: + 3ks) -3ksAk + ^^"'^ (^s - 3ks) J 
where A-^ = —{Ax + ks). The pseudoscalar masses are 



(3.13) 



^kp2 = 2 [^22 + Mn T V(^22 - Mn)2 + 471^^2], (3.14) 

where 

M22 ± Mn = -{Aj: - Sks) - SksA^ ± — -|, 

2s sm2p 

M12 = Av(As + 3A;s). (3.15) 

The pseudoscalar eigenstates are described through the mixing angle 7 given by: 

. -2Mi2 

sm27 = . (3.16) 
^(M22 - Mii)2 + 4Mf2 

For large values of s, the pseudoscalar Higgs masses are approximately given by 

ml^ -3ksA„ m% (3.17) 
^ ^ sm2/3 
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This shows that the pseudoscalar Pi is mainly Iras', whereas P2 is a mixture of ImHf 
and ImH^. Furthermore, for large values of tan/3, P2 becomes very massive. Thus Pi 
decouples from the spectrum because it is mainly a singlet and P2 decouples from the 
spectrum because it is very massive simplying that the pseudoscalar Higgs bosons would 
be impossible to produce in the limit of Yukawa unification. These qualitative features 
of the pseudoscalar spectrum survive the effects of radiative corrections. Here we find 
a qualitative distinction with the corresponding pseudoscalar spectrum of the MSSM. On 
the other hand the tree level the squared mass matrix, M|, for the neutral scalar Higgs 
bosons can be written as 

^ m|cos2/? + Asstan/3 -XsA^ + X^v^ sin 2(3 - "'^ ^ Xv2{2Xs cot P + ks ~ A^) ^ 

-XsA^ + X^v^sm2P- ^^^^^^ m| sin^ /? + ^ss cot /3 Xvi{2Xs tan f3 + ks ~ Aj:) 

y Xv2{2Xscotp + ks~ A^) Xvi{2Xstan(3 + ks- A^) (4fc^5^ + fcsQfc) - ^-4.^^^^sin2/3 j 

(3.18) 

We shall label the CP-even Higgs eigenstates of eq.(3.18) as 5*1, S2, and 5*3, respectively 
in order of increasing mass. It is not very illuminating to present complicated analytical 
expressions for them, and so we will defer the numerical results to Sec. 4. It turns 
out that in the limit of large tan/?, the lightest CP-even Higgs boson is almost a pure 
Re(i72), whereas the heaviest CP-even Higgs boson is predominantly the singlet, KeS. 
Furthermore, the lightest Higgs boson mass has an upper bound which is close to the 
upper bound in the MSSM in the same limit. This is a consequence of the fact that the 
upper bound on the lightest Higgs mass in the two models, including radiative corrections, 
differs by a term that is proportional to sin 2(3, and is, therefore, small . 



The neutralino mass matrix (M^o) in the Non Minimal Supersymmetric Standard 
model can be written as [^] : 

^ Ml —niz cos (3 sm9m mz sin (3 sin 9 w ^ 

M2 mz cos (3 cos 6w — mz sin P cos 9^^ 

—mz cos (3 sin9w mz cos P cos 9w As Aw sin /3 , (3.19) 

mz sin (3 sin 9w —mz sin (] cos 9u! As Xvcos(3 

Aw sin /3 Xvcos(3 2ks 
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where Mi and M2 are the masses of the U{1)y and SU{2)l gauginos (A'; A", a = 1, 2, 3), 
respectively, and where we have chosen the basis (— iA', — iA^, ^//i, ^_f/|) ^s)- The grand 
unification condition leads to the mass relation Mi = (5/3)M2 tan^ The gluino mass 
is M3; the gluino does not mix with the rest of the neutralinos. Because neutralinos are 
Majorana particles, the mass matrix (3.19) is in general comples symmetric, and hence 
can be diagonahzed by only one unitary matrix N, i.e., N*M^oN~^ — M^. We shall label 

the neutralino masses in the ascending order (of magnitude) as M^o, i — 1, ,5. 

In the limit of large s, we have the the approximate expressions for the (magnitude 
of) neutralino masses in increasing order: 

Ml, M2, Xs,Xs,2ks, (3.20) 

with the lightest neutralino being a predominantly a bino, which satisfies the simple mass 
relation 

Ml ~ ^^ii^Mi/2 ~ O.45M1/2. (3.21) 

All the neutralinos with the exception of the heaviest one have negligible singlet com- 
ponent; the singlet decouples from the remainder of the spectrum. We note that two of 
the neutralinos are nearly degenerate, thus leading to a pseudo-Dirac neutralino, as in 
MSSM, in the hmit of large tan /3. 

The chargino mass matrix M^± is given by: 



^ M2 V^mw sin (3 ^ 



(3.22) 



y -\/2miyCos/3 —As J 
and is same as in MSSM with the role of played by As. The masses and the composition 
of the charginos are obtained by diagonalizing the matrix (3.22) via the biunitary trans- 
formation UMy^±V~^ = M^±, where U and V are 2x2 matrices which diagonalize the 
hermetian matrices M-^±M^± and M^±M^±, respectively. Similarly, the masses of other 
sparticles are same as in MSSM with n replaced by As. 
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4 Numerical Results and Discussion 



Having described the NMSSM with large tan (3 in detail in the previous section, we now 
turn to obtaining the particle spectrum of the model numerically. 

Writing down the coupled set of the RG equations for the 24 parameters of the 



model p2| and including the contributions of hi, and h^- [|13l, the parameter space 
of the model is scanned by taking values of the input parameters (M1/2, mo, A, h, A, k) at 
Mx which are evolved down to low energies Qo to obtain the values of of the parameters 
tan/5, r = s/v, Ax, A^, ks . We also evolve the values of the soft masses appearing on 
the left hand sides of eqs. ( |3.1| ) - ( |3.3| ), and compare their values to the combinations 
of the parameters appearing on the right hand sides of these equations as obtained from 
the RG evolution. The input parameters are chosen so as to satisfy the constraints de- 
scribed in the previous section, namely the absence of electric charge breaking vacua, the 
charged Higgs mass squared remaining positive, and the SU (2) x f/(l) breaking minimum 
being energetically favorable. It turns out that the first of the minimization conditions, 
eq. (|3.1|) , is the one which is most sensitive to the choice of initial conditions reflecting 
the fine tuning discussed in Sec. 3. The parameters are chosen so as to study ri, r2 and 
r3, which are defined as the difference between the left and right hand sides of the three 
minimization eqs. (3.1) - (3.2) divided by the right hand side of each of these equations, 
and study the change in sign that these suffer as the parameters are varied. There are 
enormous difficulties in trying to achieve a simultaneous solution to = 0, i = 1, 2, 3. 
In particular = requires the presence of values for |yl|/mo of almost 3 or more. This 
requirement has a profound impact on the particle spectrum of the model. In particular, 
r, the ratio of the singlet to the doublet vacuum expectation value persistently remains 
large for the choice of parameters considered with large tan/5, corresponding to the VEV 
of s being of the order of 10 TeV. This is substantially different from the situation when 
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tan /3 ~ 1 [|T^ . We note that the singlet vacuum expectation value is not constrained by 
the experimental data. 

In Fig. 1, we show a typical evolution of the three soft SUSY breaking mass parameters 
m|^^, ^'^d m| from Mx down to the low scale Qq with a choice of parameters such 
that all constraints are satisfied, and we are in the neighbourhood of an SU{2) x f/(l) 
breaking vacuum. We note that because of the possibility of large value of mt(m() = 181 
GeV, we have a large value for h so that the Yukawa couplings dominate over the gauge 
couplings in the evolution of these parameters. This in turn forces the mass parameters 
to remain large at large momentum scales compared to their values at smaller momentum 
scales. 

In supersymmetric theories with R-parity conservation, the lightest supersymmetric 
particle generally turns out to be the lightest neutralino. In the NMSSM the neutralino 
mass matrix is given by eq.( p.l9|) and its general properties are discussed in The 
parameters that determine the mass matrix are X, k, s,tan /3, Mi and M2. Choosing the 
input parameters at Mx so that all the constraints of the Sec. 3 are satisfied, we obtain 
values for these parameters at Qq and then the neutralino mass matrix may be evaluated 
numerically. The chargino mass matrix eq. (|3.22| ) may also be evaluated in a similar 
manner. 

One result of this procedure is shown in Fig. 2, where we plot the lightest neutralino 
and chargino masses for a specific choice of input parameters of Table 1, as obtained from 
RG evolution, as a function of the top quark mass. We have found from our scan of the 
parameter space, that the lightest neutralino is almost a pure bino in the limit of large 
tan jS. Furthermore, all other neutralinos except the heaviest one have a negligible singlet 
component, indicating that the singlet completely decouples from the lighter neutralino 
spectrum. These properties of the neutralino spectrum are shown in the second and third 
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columns of Table 2. The mass of the lightest neutralino in this case is determined by 
the simple mass relation for the bino Mi ~ ^1/2 — 0.45Mi/2- The masses of the 

heavier neutralinos lie in the range of 0.5 — 1 TeV. Furthermore, the lightest chargino mass 
bears a relation to M1/2 similar to the neutralino relation, with ai in eq. (3.21) replaced 
by 0^2, reflecting that it is primarily a charged wino. This is seen from the fourth, fifth 
and sixth columns of Table 2. The heavier chargino mass is found to be ~ 1 TeV. We 
further note that two of the neutralinos are nearly degenerate, and lead to a pseudo-Dirac 
neutralino. Also the second lightest neutralino is primarily a wino and degenerate in mass 
with the lightest chargino. These characteristics are similar to those found in MSSM. The 
gluino mass is found to be 1.6 TeV for the choice of parameters of Fig. 2 and follows from 
a relation similar to eq. (12) with ai replaced by as- 

We now come to the spectrum of CP-even Higgs bosons of the model. In order 
to understand the quantitative features of the results we have obtained for the lightest 
CP-even Higgs boson, we need to go into some detail regarding the actual choices of 
parameters entering the computation and the correlations between the various elements 
of the spectrum. We have already reviewed the corresponding situation in the case of the 
MSSM in Sec. 2, and here we will discuss the comparison of the NMSSM with that of the 
MSSM. In Fig. 3, we plot for typical and reasonable values of the input parameters, in the 
region where the vacuum is expected to lie, the mass of the lightest CP-even Higgs boson 
as a function of the top-quark mass mt{mt) in the range that is most favoured under these 
boundary conditions [0, ^ 0, |^, Q . The choice of parameters here is closely related to the 
family of solutions studied extensively in Ref. [0, and would serve as a typical example 
of the numbers we have explored. In the MSSM when the mass of its unique CP-odd 
Higgs boson ^ m^, the substantive part of the radiative correction is picked up by 
the lighter of the CP-even bosons, As tua approaches mz, the radiative corrections 
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are shifted to the heavier of the CP-even Higgs bosons, H^. Such a feature is observed 
here as weU: for those choices of parameters in Table 1 that yield a somewhat smaller 
(the mass of lightest CP-odd Higgs boson in NMSSM), we find that the radiative 
corrections to the lighest CP-even Higgs, 5*1, are smaller. We note, however, that the 
CP-odd Higgs bosons here are always massive in the limit of large tan (3 in contrast to the 
situation that prevails in the MSSM, where tua even in the vicinity of mz is plausible |0. 
Due to the complexity of the system under investigation and the difficulty in controlling 
the numerical choice of the parameter A for a given h, with all other parameters held 
fixed, we do not know how precisely close the choice of parameters of Table 1 are to a 
genuine ground state. Furthermore, we note that the clarity with which the correlations 
have been observed between rriA and rrifiO in MSSM do not have a simple parallel here due 
to the presence of a larger number of physical states. A more precise, albeit prohibitively 
time consuming, determination could then ensure that the spurious wobble seen in Fig. 
3 is eliminated, and would establish a more reliable correlation between increasing h and 
the rise of the mass of 5*1 and the correlations with mp^. Furthermore, a refinement of 
the choice of parameters, based on the minimization of the one-loop effective potential, 
could stablize the figures presented here. 

We note that the lightest Higgs bosons mass ~ 130 GeV for a wide range of parame- 



ters which nearly saturates the upper bound of 140 GeV |JT3[, and lies in the same range 
as in MSSM with large tan/3. This is a consequence of the largeness of tan/5: the con- 
tribution to the tree level mass which depends on the trilinear coulings A is small, being 
proportional to sin^ 2(3, so that the upper bound on the lightest Higgs mass reduces to the 
corresponding upper bound in MSSM when the appropriate identification of the param- 
eters is performed 0]. We also note that the upper bound on the lightest Higgs mass 
depends only logarithmically on r, and hence on the singlet vacuum expectation value s, 
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in the limit of large r, which, therefore decouples from the bound [^. Furthermore, the 
lightest Higgs boson in almost a pure doublet Higgs (Re -^2); "with the singlet component 
being less that 1% in the entire range of parameters considered. It is only the second 
heavier CP even Higgs boson S2 that is predominantly a singlet. Its mass ranges between 
740 GeV and 2.3 TeV. The heaviest CP even Higgs boson 6*3 is again predominantly a 
doublet Higgs (Re H^) with its mass varying between 4 — 6 TeV. This implies that all 
the CP-even Higgs bosons, except the lightest one, decouple from the spectrum. These 
features of the spectrum of the CP-even Higgs bosons of NMSSM for large values of tan /3 
are clearly seen from Table 3, where we show the mass and composition of Si for a wide 
range of input parameters. The results presented above, that the lightest Higgs boson is 
almost purely a doublet Higgs at large tan/5, are in contrast to the situation with low 
values of tan (3, where the lightest CP-even Higgs boson contains a large admixture of the 
gauge singlet field S [|12[ |2^ . 

In Fig. 4 we plot for the typical values of the input parameters of Table 1 the mass 
of the lightest CP-odd Higgs boson as a function of mt(mj). We note from Table 3 that 
both CP-odd Higgs bosons Pi and P2 are heavy, their masses being in the range of 2 TeV 
and 6 TeV, respectively. Also, the lightest CP-odd state is predominatly a Higgs singlet, 
thereby effectively decoupling from the rest of the spectrum. The charged Higgs boson 
mass rric lies, for most of the cases that we have studied, in the range 1 — 2 TeV. 

In order to discuss the features of the the sfermion spectrum, we first recall some of 
the features of the spectrum of the MSSM. In the MSSM it has been observed ||^ that 
the presence of large Yukawa couplings for the b-quark as well as the r lepton, as well as 
the presence of large trilinear couplings, could lead to the lighter of the scalar r's tending 
to become lighter than the lightest neutralino, which is the most favoured candidate in 
such models for the lightest supersymmetric particle. In particular, in order to overcome 
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cosmological constraints for given values of M1/2, lower bounds on mo were found to 
emerge. In turn, increasing ttiq implies ever decreasing (clear correlations have been 
described for the case of A = in Ref. [^j) thus leading to further constraints on the pa- 
rameter space of the MSSM. For large values of tan (3 the competing tendencies between 
the lighter scalar tau mass and rriA play an important role in MSSM in establishing a 
lower bound ~ 450 GeV on M1/2. Given the complexity of the system of equations, it has 
not been possible to to extract similar lower bounds on M1/2 in NMSSM. Nevertheless, in 
Ref. [jl3| the intimate link between the ground states of the two models has been estab- 
lished and a much more sophisticated and time consuming analysis of the present model 
is also likely to yield a lower bound that is unlikely to be very different from the one 
obtained in MSSM. As a result, in confining ourselves to numbers of this magnitude and 
higher, we find a heavy spectrum. More recently |^ further experimental constraints on 
MSSM have been taken into account resulting in an extension of the minimal assumptions 
at Mx by including non-universality for scalar masses. Indeed, in the present analysis 
similar problems have been encountered with some of the choice of parameters studied 
in Ref. [|13|, with rrif^ tending to lie below the mass of the lightest neutralino due to the 
persistent presence of large Yukawa couplings and more so due to the large trilinear cou- 
plings dictated by eq.( p.3|) . Nevertheless, given the fact that the present work minimizes 
the tree-level potential and that the violations of cosmological constraints are not seri- 
ous, in that minor adjustments of |A|/mo solve this problem efficiently, we consider the 
regions of the parameter space we have explored to be reasonable. Furthermore, it could 
be that the extension of minimal boundary conditions along the lines of Ref. could 
provide alternative and elegant solutions to this problem, while preserving the existence 
of relatively light scalar r's as a prediction of the unification of Yukawa couplings in the 
NMSSM as well as in the MSSM. 
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The heaviest sfermions in the spectrum of NMSSM, as in MSSM, are the scalar quarks, 
which tend to be much heavier, in the TeV range. The S'O(IO) property that the scalar 
b-quarks are as massive as the scalar top-quarks is preserved in the NMSSM. 

To summarize, through a detailed analysis of NMSSM presented here, we have shown 
that all the particles, except the lightest CP - even Higgs boson, implied by supersymmetry 
are heavy for large values of tan (3. The gauge singlet field S decouples, both from the 
hghtest Higgs boson as well as the neutrahnos. The LSP of the model continues to be, as 
in MSSM, the lightest neutralino that is primarily a bino in composition, with the lighter 
scalar r having a mass in the neighbourhood of the LSP mass. The remainder of the 
spectrum tends to be heavy, from 1 to a few TeV. We note that the NMSSM in the large 
tan 13 regime rests on a delicately hinged system of equations and constraints. Although 
it provides a good testing ground for the stability of the predictions of the MSSM, in 
practice it deserves great care in its treatment. 
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Table Captions 

Table 1 Sample of values of input parameters (M1/2, mo, A, h, k, A) and the computed 

values of different parameters {■mt{mt), tan/3, r. Ax, Ak and ks). All mass parameters are 
in units of GeV. 

Table 2 Masses and compositions of neutralino x° and chargino xf states for the values 
of input parameters of Table 1. All masses are in GeV. 

Table 3 Masses and compositions of the CP-even (Si) and CP-odd (Pj) Higgs bosons. 
The composition of CP-even Higgs bosons is in terms of the basis {ReH'^, ReH2, ReS), 
whereas the basis for CP-odd Higgs bosons is (^4°, ImS). All masses are in GeV. 

Figure Captions 

Fig. 1 The evolution of soft supersymmetry breaking mass parameters from the grand 
unified scale Mx to Qq defined in the text. The input parameters are M1/2 = —700, mo 
and A = 1600 (all in GeV). The other parameters are h = 1.5, A = 0.40 and k = -0.10. 
The associated value of the top quark mass is 181 GeV. 

Fig. 2 The lightest neutralino and chargino masses as a function of mt{mt). The input 
parameters are M1/2 = —700, mo = 800, A — 2200 (all in GeV), with the remaining 
parameters varied to guarantee a solution. 

Fig. 3 The lightest CP-even Higgs boson mass as a function of mj. The range of 
parameters is as in Fig. 2. 

Fig. 4 The hghtest CP-odd Higgs boson mass as a function of rrit. The range of param- 
eters is as in Fig. 2. 
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# 


Ml/2 


mo 


A 


h 


A 


k 


mt{mt) 


tan [3 


r 




Ak 


ks 


1 


-700 


800 


2200 


0.75 


0.1 


-0.1 


170.5 


54 


85 


908 


2131 


-1453 


2 


-700 


800 


2200 


1.00 


0.2 


-0.1 


176.5 


58 


58 


807 


2100 


-975 


3 


-700 


800 


2200 


1.25 


0.3 


-0.1 


179.6 


60 


49 


745 


2070 


-807 


4 


-700 


800 


2200 


1.50 


0.4 


-0.1 


181.4 


62 


44 


703 


2045 


-717 


5 


-700 


800 


2200 


1.75 


0.4 


-0.1 


182.5 


63 


51 


698 


2067 


-835 


6 


-700 


800 


2200 


2.00 


0.5 


-0.1 


183.3 


64 


47 


675 


2048 


-755 


7 


-700 


800 


2200 


2.25 


0.6 


-0.1 


183.8 


64 


44 


657 


2032 


-700 



Table 1 
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# 












1 


-309 
-579 
877 
-887 


0.99,-0.01,-0.06,0.02,0.00 
0.02,0.98,0.16,-0.11,0.00 
0.03,-0.04,0.70,0.70,0.00 

0.05,-0.20,0.70,-0.70,0.00 


579 


0.99,0.15 


0.97,0.22 




-2906 


0.00,0.00,0.00,0.00,0.99 


886 


-0.15,0.99 


-0.22,0.97 


2 


-309 
-580 
941 
-949 


0.99,-0.01,-0.05,0.02,0.00 
0.01,0.99,0.14, -0.09,0.00 
0.02,-0.04,0.70,0.70,0.00 

0.05,-0.16,0.70,-0.70,0.01 


581 


0.99,0.12 


0.97,0.22 




-1950 


0.00,0.00,-0.01,0.00,0.99 


949 


-0.12,0.99 


-0.22,0.97 


3 


-309 
-581 
936 
-944 


0.99,-0.06, -0.05,0.02,0.00 
0.01,0.97,0.14, -0.09,0.00 
0.02,-0.04,0.70,0.70,0.00 

0.05,-0.16,0.69, -0.70,0.00 


581 


0.99,0.10 


0.98,0.19 




-1613 


0.00,0.00,-0.02,0.01,0.99 


945 


-0.10,0.99 


-0.19,0.98 


4 


-309 
-582 
1001 
-1007 


0.99,0.00,-0.05,0.02,0.00 
0.01,0.99,0.12,-0.07,0.00 
0.02,-0.03,0.70,0.70,0.00 
0.04,-0.14,0.70,-0.7,0.04 


582 


0.99,0.12 


0.99,0.07 




-1436 


0.00,0.00,-0.03,0.02,0.99 


1008 


-0.12,0.99 


-0.07,0.99 


5 


-309 
-582 
979 
-986 


0.99,0.00,-0.05,0.02,0.00 
0.01,0.99,0.13,-0.08,0.00 
0.02,-0.03,0.70,0.70,0.00 
0.05,-0.14,0.70,-0.7,0.00 


582 


0.99,0.10 


0.99,0.17 




-1669 


-0.33,0.00,-0.02,0.00,0.99 


1018 


-0.10,0.99 


-0.17,0.98 


6 


-309 
-582 
978 
-984 


0.99,0.00,-0.05,0.02,0.00 
0.01,0.99,0.13,-0.08,0.00 
0.02,-0.04, 0.70,0.70,0.00 
0.05,-0.14,0.70,-0.7,0.00 


581 


0.99,0.10 


0.94, 0.35 




-1509 


0.00,0.00,-0.02,0.01,0.99 


985 


-0.10,0.99 


-0.35,0.94 


7 


-309 
-582 
998 
-1004 


0.99,0.00,-0.05,0.01,0.00 
0.01,0.99,0.12,-0.07,0.00 
0.02,-0.03,0.70,0.70,0.00 
0.05,-0.14,0.70,-0.7,0.04 


582 


0.99,0.10 


0.98,0.17 




-1402 


0.00,0.00,-0.03,0.03,0.99 


1004 


-0.10,0.99 


-0.17,0.98 



Table 2 



# 




Com|)<)sitiou 


irir, 


Com|")<)sition 


1 


131 


-0.02,0.99,0.00 


3048 


0.00,1.00 




2314 


0.00,0.00,0.99 








6039 


0.99,-0.02,0.00 


6918 


1.00,0.00 


2 


132 


0.02,0.99,-0.02 


2479 


0.00, 1.00 




1326 


0.00,0.02,0.99 








4822 


0.99,-0.02,0.00 


6216 


1.00,0.00 


3 


130 


0.02,0.99,-0.05 


2238 


0.00, 1.00 




967 


0.00,0.05,0.99 








4342 


0.99,-0.17,0.00 


5990 


1.00,0.00 


4 


119 


0.02,0.99,-0.09 


2099 


0.00, 1.00 




774 


0.00,-0.9,0.99 








4297 


0.99,-0.02,0.00 


6394 


1.00,0.00 


5 


133 


0.02,0.99,-0.04 


2275 


0.00, 1.00 




1031 


0.00,-0.04, 0.99 








5116 


0.99,-0.02,0.00 


6800 


1.00,0.00 


6 


128 


0.02,0.99,-0.07 


2153 


0.00,1.00 




857 


0.00,0.07,0.99 








4823 


0.99,-0.02,0.00 


6636 


1.00,0.00 


7 


119 


0.02,0.99,-0.01 


2067 


0.00, 1.00 




739 


0.00,0.10,0.99 








4644 


0.99,-0.02,0.00 


6557 


1.00,0.00 
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